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Abstract 



We consider supersymmetric inflation models in which inflation occurs at an inter- 
mediate scale and which provide a solution to the \i problem and the strong CP 
problem. Such models are particularly attractive since inflation, baryogenesis and the 
relic abundance of cold dark matter are all related by a set of parameters which also 
affect particle physics collider phenomena, neutrino masses and the strong CP prob- 
lem. For such models the natural situation is a universe containing matter composed 
of baryons, massive neutrinos, lightest superpartner cold dark matter, and axions. 
The present day relic abundances of these different forms of matter are (in principle) 
calculable from the supersymmetric inflation model together with a measurement of 
the CMB temperature and the Hubble constant. From these relic abundances one 
can deduce the amount of the present day dark energy density. 
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1 Introduction 



Recent data on the cosmic microwave background (CMB) radiation |J provides strong 
support for inflation, by measuring the first, second and third peaks of the angular 
power spectrum. Inflation therefore seems to be increasingly well established. Low 
energy supersymmetry (SUSY) is perhaps less well established (though there is con- 
siderable indirect evidence for it [0]), but is certainly desirable from many points of 
view, and has the advantage when combined with inflation of helping to ensure that 
the inflaton potential is sufficiently and naturally flat j3|. The CMB data also sup- 
ports a ACDM universe in which the energy density is dominated by dark energy 
(DE) (corresponding possibly to a cosmological constant A), and cold dark matter 
(CDM) §. ' 

In this paper we consider supersymmetric inflation models in which inflation oc- 
curs at an intermediate scale || [7], [|, and which provide an intermediate scale 
solution to the \i problem and a solution to the strong CP problem via the Peccei- 
Quinn (PQ) mechanism ||. When right-handed neutrinos are included, such models 
may give baryogenesis via leptogenesis. They will have calculable CDM and DE relic 
densities. And they may have supergravity (SUGRA) mediated SUSY breaking with 
no moduli or gravitino problems. A main point of the paper is to show that in such 
models in which all these features are simultaneously present there will be fewer pa- 
rameters than in models in which these problems are separately addressed. We will 
discuss an explicit example of such a model in order to demonstrate the connections 
between the physics of each of the separate sectors, and the resulting enhanced pre- 
dict ivity leading for example to connections between collider physics experiments and 
cosmological observations. 

Within such a framework one may expect on general grounds that our Universe 
contains sizeable relic abundances of baryons (from leptogenesis), axions (a) (from 
the solution to the strong CP problem), as well as weakly interacting massive parti- 
cles (WIMPS). In R-parity conserving SUSY the WIMP is identified as the lightest 
supersymmetric particle (LSP), and this is often assumed to be the lightest neu- 
tralino of the minimal supersymmetric standard model (MSSM). We shall argue from 
this perspective that the LSP could equally well be a lighter stable singlet (singlino) 



||T0| identified with an axino |TT| or inflatino ||12|| . Neutrino masses provide smaller 
amounts of hot dark matter. Within this approach all these forms of matter in the 
universe will have calculable relic abundances, given a measurement of the CMB 
temperature and the Hubble constant, which are related to the parameters of the 
underlying supersymmetric theory which may be determined from particle physics 
experiments. This amounts to a generalisation of the observation made a few years 
ago that the relic abundance of neutralinos is related to the parameters of the super- 
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symmetric theory. Further, a given parameter typically is relevant to more than one 
relic abundance, so the total number of parameters is fewer than when the various 
forms are considered separately. From the calculated present day relic abundances 
of matter one can deduce the amount of the present day DE density, even without 
specifying the physics of DE. 



2 The Omega Problem in Supersymmetric Infla- 
tion 

In general the ratio of the total density of the universe ptot to critical density p CT u is 
given by Vt to t where 

^tot = ^7 ^matter 

+ tt DE (1) 

and f2 7 , flmatter, &de are the ratios of radiation density p 7 , matter density p ma tter 
and DE density pde to critical density p cr it, and the radiation density is unimportant 
P-y "C Pmatter- Note that the critical density is a function of time and in the present 
epoch p crit = 3M 2 P Hl = (3/i 1 / 2 x 1(T 3 eV) 4 ~ (M^/M P ) 4 where M w is the weak 
scale, Mp is the Planck scale, and H = 100/ikm.s _1 Mpc _1 is the present day Hubble 
constant, with h = 0.7 ± 10% |13[. From observation || VLde ~ 2/3 while flmatter ~ 



1/3 and VL to t is very close to unity. Inflation predicts Vt to t = 1, for all times after 
inflation. The matter contributions consist of (at least) 

^matter = + &d + ^LSP + ^a (2) 

The most recent data Q is consistent with nucleosynthesis estimates of fi& ~ 0.04, 
where the baryons (b) in the universe are mainly to be found in dark objects. The 
determined value of SIqdm ~ 0.3 contains unknown relative contributions from CIlsp 
and Q a . Super-Kamiokande sets a lower limit on neutrino masses Z)i m ^ > 0.05 eV 
which corresponds to Q u > 0.003. In hierarchical neutrino mass models the lower 
bound is saturated and interestingly the neutrino density is then comparable to the 
visible baryon density fl s tars ~ 0.005. 

We now discuss how to calculate the relic densities from our present rather prim- 
itive state of knowledge about the comprehensive theory needed to really do that. [] 
In the following it is important to keep in mind that our present inability to calculate 
the ratios of different forms of matter should be distinguished from the ability to 
calculate them in principle. Later we discuss a simple model which illustrates these 
ideas. 



lr rhe calculation will involve non-pcrturbative cosmological effects during the reheating process, 
such as preheating @ |l|, [H| . 
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Of the three likely dominant densities f£&, £Ilsp arid Q a , the latter arises from 
non-relativistic axions being produced at the QCD scale by the usual misalignment 
mechanism, and is the most difficult of the three to estimate since it depends on a 
randomly selected angle. Q Therefore we shall focus mainly on the question of how 
to calculate flj, and Qlsp- 

The present day value of Qb corresponds to the ratio of baryon number density 
to entropy density of the universe Yj, = n^/s ~ 0.7 x 1CU 10 , assuming ng = 0. For 
the calculation of Y&, we shall restrict ourselves here to the so-called leptogenesis 
mechanism. The basic idea of leptogenesis Jl8|, [19|] is that right-handed neutrinos (or 
possibly sneutrinos) are copiously produced in the early universe, then decay to pro- 
duce lepton number (and hence B — L) asymmetry. The lepton number asymmetry 
is subsequently converted into baryon number asymmetry by sphaleron interactions. 
Concerning the three Sakharov conditions: the CP-violation originates from complex 
Yukawa coupling constants fj lepton number violation originates from the Majorana 
mass of the right-handed neutrinos, and baryon number violation from sphalerons. 
Concerning the out-of-equilibrium condition, in the conventional approach to leptoge- 
nesis it is assumed that the right-handed neutrinos are produced by their couplings to 
other particles in the thermal bath, but that these couplings are sufficiently weak that 
the decays occur out-of-equilibrium, leading to a narrow range of couplings JJO, EI |. 



From the perspective of inflation the conventional leptogenesis picture will change if 
the reheat temperature is below the mass of the lightest right-handed neutrino. In 
this case right-handed neutrinos may be produced during reheating via (direct or in- 
direct) couplings to the inflaton field, and may then be produced with masses greatly 
exceeding the reheat temperature, providing only that they are lighter than the in- 
flaton field. In this case the out-of-equilibrium condition is automatically satisfied 
during reheating. This second mechanism is preferred from the point of view of the 
gravitino constraint, since in this case if the reheat temperature is below the limit 
Tr < 10 9 GeV then thermally produced gravitinos are not a problem.^ 

In the conventional case the LSP is regarded as the lightest neutralino xi ( a 
linear combination of neutral bino B, wino W3, and higgsinos H±, H2) and it is pro- 
duced in thermal equilibrium at temperatures above its mass ~ Mw- According to 
the following very crude argument when their annihilation rate r ~ arif ~ rif/M^ 
becomes smaller than the expansion rate of the universe H ~ Tj/Mp (at a tem- 
perature Tf ~ Mw), they freeze out of the thermal bath, and the present day 
matter density is then p x ~ M w nf(T/M w ) 3 ~ M W T 3 /M P including the dilution 

2 Any relativistic axions produced during (p)reheating will be red-shifted away. 
3 Complex SUSY mass parameters are not relevant for leptogenesis if its scale is above that of 
SUSY breaking. 

4 If the gravitino is the LSP then this limit on the reheat temperature may be relaxed |E3 . 
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factor (Rf/R) 3 ~ (T/Mw) 3 - The current temperature is obtained by roughly equat- 
ing p x ~ p 7 m T 4 which gives T ~ M^/Mp. Inserting this temperature we find 
p x ~ (M^r/Mp) 4 ~ Pent- Many more careful analyses have been performed in order 
to obtain precise estimates for p x by considering the detailed annihilation channels 
within different regions of SUSY parameter space. What concerns us more here is 
how including the effects of inflation will change this simple framework. As in the 
discussion of leptogenesis (above) from the perspective of inflation and baryogenesis, 
and in general a broader picture, the conventional picture will change if the reheat 
temperature is below the mass of the lightest neutralino, and they are produced by 
non-thermal processes during the reheating process after inflation. Another way in 
which the physics might differ is if the LSP is not the lightest neutralino, but instead 
some lighter singlino associated with an axino or inflatino. For example, if the lightest 
neutralino freezes out then decays into an axino (a) then the present day axino density 
is suppressed by the ratio of axino mass to lightest neutralino mass p a = (m a /m x )p x 
In order to learn the actual relic density of LSPs (or axions or any candidate) it 



is necessary to actually calculate it; detecting LSPs or axions is possible even if the 
relic density is well below fl ma tter of order 1/3 p3| . 



3 An Example 

3.1 Why Intermediate Scale Inflation? 

We now wish to consider a specific example of a model which addresses the particle 
physics issues mentioned earlier, in order to illustrate many of the general features 
that we have discussed above. The brand of inflation most closely related to particle 
physics seems to be hybrid inflation which may occur at a scale well below the Planck 
scale, and hence be in the realm of particle physics p4[ . The next question is what 



is the relevant scale at which hybrid inflation takes place? One obvious possibility is 
to associate the scale of inflation with some grand unified theory (GUT) symmetry 



breaking scale, as originally conceived by Guth p5 ]. However it is somewhat ironic 
that hybrid inflation at the GUT scale faces the magnetic monopole problem, which 
was precisely one of the original motivations for considering inflation in the first place! 
Although in certain cases this problem may be resolved, there are typically further 
symmetry breaking scales below the GUT scale at which discrete symmetries are 
broken, leading to problems with cosmological domain walls. 

As an example of the difficulties faced by GUT scale inflation models, consider 
the breaking of the Pati-Salam symmetry group SU(4) x SU(2) x SU(2) down to 
the standard model gauge group. The minimal symmetry breaking potential in this 
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model 



26[ contains a singlet which could be a candidate for the inflaton of hybrid 
However it was immediately realised that such a scenario would face a 



inflation. 

magnetic monopole problem since the gauge group is unbroken during inflation, and 
only broken by the choice of vacuum after inflation [26|. A possible solution to this 
problem is to consider the effect of higher dimension operators in the superpotential 
which for a range of parameters have the effect of breaking the gauge symmetry 
during inflation |27|]. Such a scheme can also address the /i problem and the strong 



CP problem, at the expense of introducing additional singlets which develop vacuum 
expectation values (VEVs) at the PQ symmetry breaking scale |§. However in this 
scheme it turns out that the vacuum in which PQ symmetry is unbroken is preferred 
below the reheat temperature of the model so that PQ symmetry is not broken after 
inflation ||27|| . Moreover PQ symmetry breaking is associated with the breaking of 



discrete symmetries in the model, so that it would lead to the domain wall problem 
in any case |27|. The only solution to these problems seems to be to assume that PQ 
symmetry is also broken during inflation, but since the inflaton has zero PQ charge, 
and since the inflation scale of order 10 14 GeV [27| and hence much larger than the 
PQ symmetry breaking scale, this assumption seems rather questionable. 

Intermediate scale hybrid inflation immediately solves both the magnetic monopole 
problem and the domain wall problem. The idea is simply that there is a period of 
hybrid inflation occuring below the GUT scale at the PQ symmetry breaking scale 
itself, in which the inflaton carries PQ charge and so the choice of domain is fixed 
during inflation. The universe therefore inflates inside a particular domain, and the 
magnetic monople relics produced by the GUT scale symmetry breaking are inflated 
away. This provides a powerful motivation for intermediate scale inflation, which is 
is the subject of this paper. We now turn to an explicit example of an intermediate 
scale inflation model. 



3.2 Intermediate Scale Supersymmetric Inflation Model 

The model we consider || is a variant of the NMSSM. This model has a SUGRA 
foundation ||, and leptogenesis and reheating has been studied [0, and preheating 



17| , |2q] has been demonstrated not to lead to over-production of either axions or 
gravitinos. The model provides a solution to the strong CP problem and the \i 
problem, with inflation directly solving the monopole and domain wall problems at 
the inflation scale. It is therefore a well motivated, successful model that has been 
well studied and does not appear to suffer from any embarrassing problems, and is 
therefore a suitable laboratory for our discussion here. This variant of the NMSSM 
has the following superpotential terms involving the standard Higgs doublets H u , Hd 
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and two gauge singlet fields (inflaton) and N, 

W = XNH u H d - k(f)N 2 (3) 

where A, k are dimensionless coupling constants. Notice that the standard NMSSM is 
recovered if we replace the inflaton by N. However this leads to the familiar domain 
wall problems arising from the discrete Z 3 symmetry. In this new variant, the Z 3 
becomes a global Peccei-Quinn U(1)pq symmetry that is commonly invoked to solve 
the strong CP problem ||. This symmetry is broken in the true vacuum by inter- 
mediate scale and N VEVs, where the axion is the pseudo-Goldstone boson from 
the spontaneous symmetry breaking and constrains the size of the VEVs. With such 
large VEVs this model should be regarded as giving an intermediate scale solution to 
the \i problem, and as such will have the collider signatures discussed in | fTU| . 

We can make the 0-field real by a choice of the (approximately) massless axion 
field. We will now regard and N to be the real components of the complex singlets 
in what follows. When we include soft SUSY breaking mass terms, trilinear terms 
Akk<f)N 2 + h.c. (for real Ak) and neglect the \NH u Hd superpotential term, we have 
the following potential: 

V = V + k 2 N* + -m 2 ((f>)N 2 + -ml<f> 2 (4) 

where m 2 (0) = m? N + 4k 2 <p 2 — 2kA k (j). We can identify the various elements of the 
potential: Vq arises from some other sector of the theory, SUGRA for example, and 
dominates the potential ||; the soft SUSY breaking parameters A k and are 
generated through some gravity- mediated mechanism with a generic value of O {TeV) . 

The basic idea of hybrid inflation is very simple ||. For the field dependent mass 
of the N scalar positive, m 2 (0) > 0, then N = since its potential has positive 
curvature. With N = the potential becomes very simple indeed, 

V = V + l -ml4> 2 (5) 

We shall assume that comes from no-scale SUGRA, and vanishes at the Planck 
scale ||, so that it is generated through radiative corrections such that m 2 ~ —k 2 A\. 
Since m 2 is negative, during inflation is slowly rolling away from the origin, and 
therefore we have inverted hybrid inflation. When exceeds a critical value 

C = (A k /Ak)(l - sfl - Am 2 N /A 2 ) (6) 

the sign of the field dependent mass of iV will become negative, m 2 (0) < 0, and the 
iV field will no longer be pinned to zero, but will roll out to the global minimum of 
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the potential, 



(<j ) ) = A k /Ak, (N) = (A k /2y/2k)yJl-4:m 2 N /Al (7) 

Together with our assumptions about the SUSY parameters, this implies 

A k ~ k(f) c ~ k(N) ~ k((/>) ~ X(N) = 11 ~ 10 3 Gel/ (8) 

During inflation the inflaton field must satisfy the slow roll conditions e,r] <C 1. R 
The COBE normalisation 5h = 1-95 x 10~ 5 requires the value of the inflaton mass 
from radiative corrections m| ~ —k 2 A\ ~ — (lOOeU) 2 and from Eq.|8] this implies 
that A, k ~ 10~ 10 and (AT) ~ ((f)) ~ 10 13 GeV. We address the smallness of A, k in the 
next section. The spectral index n which relates to the power spectrum P k oc A;™ -1 , 
is given by \n — 1| = 2i] — 6e ~ 2i] ~ 10~ 12 which provides a basic prediction of the 
model. The present value of the spectral index in the range 0.80 < n < 1.05 at 68% 
C.L. slightly disfavours the accurate prediction that n = 1.00 but only at the la level, 
and we may have to wait for the results from the Planck satellite which will measure 
it to an accuracy of An = ±0.01 to definitively test this prediction |29 . 



It is non-trivial that a set of parameters exists that is consistent with axion and 
SUSY physics and allows the correct COBE perturbations to be achieved by radiative 
corrections. Without SUSY one would be free to add soft scalar masses at will, but 
with SUSY one must rely on the theory which either generates soft masses of order 
a TeV, or sets them equal to zero as in no-scale SUGRA, in which case the radiative 
corrections, which are under control in the case of SUSY predict the relevant value 
of the soft parameters, without any further adjustable parameters. Thus SUSY is 
playing a crucial role in the model which is why we refer to it as a Supersymmetric 
Inflation Model. 



3.3 Role of Non-renormalisable Operators 

The couplings A, k should be thought of as effective couplings arising from non- 
renormalizable operators || |j so they are actually couplings of order unity times 
small factors arising from ratios of VEVs to the Planck mass to some power, as can 
occur in models, and are not unnaturally small. In the original model it was suggested 
that the superpotential in Eq.[| arose as an effective theory from a non-renormalisable 

5 Recall that e = \M%{V /Vf with \r)\ = \M*V"/V\, where V'(V") are the first (second) 
derivatives of the potential and Mp = Mp/8w is the reduced Planck mass. The COBE normalisation 
is given by S 2 H = (l/1507r 2 )Vb/(M^»e). 
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superpotential which at leading order is given by 

MM ,. 9 M 2 (MM) 3 ,(iV¥) 5 (M¥) 2 

= - + «hsr + d w + " ' (9) 

where two extra singlets M and M have been introduced which develop VEVs by 
a radiative mechanism < M >=< M >~ 10 14 GeV, as a result of which we recover 
the original superpotential in Eq|| and we reinterpret our couplings A, k as effective 
couplings given by 



XuX <M^>^ ^li^, 10 -., (10 ) 



Thus the underlying coupling constants A', k' are of order unity, although for the 
most part we find it convenient to discuss the model in terms of effective couplings 
A, k. The underlying theory respects a Z 3 x Z 5 symmetry, and the U(1)pq symmetry 
arises as an approximate effective symmetry, leading to an explicit contribution to the 
axion mass from the term proportional to d which tilts the axion potential slightly, 
and perturbs the 8 angle by an amount A6 w 10~ n , thereby preserving the PQ 
solution to the strong CP problem, but providing a prediction for the next generation 
of electric dipole moment (EDM) experiments. 



3.4 The Cosmological Constant Problem 

Notice that the SUGRA-derived potential contribution Vo exactly cancels with the 
other terms (by tuning) to provide agreement with the observed small cosmologi- 
cal constant. Thus we assume that at the global minimum V((4>), (N)) = which 
implies that Vq = k 2 (N) A . The height of the potential during inflation is therefore 

Since the approach has a consistent way to set the large 
cosmological constant to zero, the absence of a real solution to this problem may not 
be an obstacle to the implications of the approach. 



3.5 Parameter Counting and Singlino Mixing 

A relevant parameter count at this stage reveals two superpotential effective param- 
eters (A and k), the two soft SUSY breaking parameters (Af- and mjv), plus the 
constant energy density Vq. From these five parameters we have inflated the universe 
with the correct COBE perturbations, provided a \i term of the correct order of mag- 
nitude and solved the strong CP problem. They also govern the phenomenology of 



S 



the singlet Higgs and Higgsino components of and N which may weakly mix with 
the MSSM superfields H u , H d . For example the Higgsino mixing matrix in the basis 
H d , H u , N, 4> is 



( -A < iV > -A < # u > ^ 

-A < N > -A < H d > 

-A < H u > -A < H d > 2k<<p> 2k< N > 

V 2k <N > ) 



(11) 



The LSP will be the lightest eigenvalue of the full "ino" matrix, extended in the usual 
way to include gaugino-higgsino mixing. Clearly if k < A/2 then a singlino will be 
the LSP. In our case the singlino may be regarded as a linear combination of axino 
and inflatino. The coupling of the superfield 5* containing the singlino S is given from 
the usual result |L0] W = /i(l + eS/v)H u H d with /z = A < N >. Here e ~ v/ f a 



where the f a is the axion decay constant and v is an electroweak VEV. Thus we have 
f a ~< iV > so e ~ A, and hence 

W~\<N>(1 + \S/v)H u H d . (12) 

As usual in models based on an intermediate scale solution to the \i problem [EL [| [TIJ 
the coupling of the singlino to the neutralinos means that S nearly decouples. However 
the conservation of R-parity means that eventually the lightest neutralino produced 
in colliders must decay into the singlino, and all the collider signatures discussed in 
|T0|1 may apply. In the case that the lightest neutralino leaves the detector before it 



decays into the singlino, there will be no unconventional collider signature. In this 
case the knowledge concerning a lighter singlino will come from cosmology since the 
LSP relic density gets diluted by the ratio of the singlino to lightest neutralino masses, 
and direct dark matter searches will not see anything since the singlino LSP will not 
scatter off nuclei. 

One of the main things we want to emphasize is the connection between the 
calculation of relic densities and the other physics, via their common parameters. 
This is summarised in Table 1 for the particular model we have been discussing. The 
same parameters that control the ino mass matrix will also be involved in reheating of 
the universe after inflation, and giving the relic densities of LSP and in leptogenesis 
as we discuss in the next section. Different models may have different mechanisms to 
solve some of the problems, different reheating and preheating, and so on, but will 
still lead to a version of Table 1. 
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Table 1: This table illustrates the fact that a particular parameter of the model (columns) simulta- 
neously controls several different aspects of particle physics and cosmology (rows) which are thereby 
related. C so ft contains Ak,Tn 2 N and the other soft parameters, Ly u k contains the Yukawa coupling 
constants controlling all fermion masses and mixing angles. 



3.6 Preheating /Reheating 

It is usually assumed that inflation ends with the singlets <f), N oscillating about their 
global minimum. Although the final reheating temperature is estimated to be of order 
1 GeV [P], during the reheating process the effective temperature of the universe, 
as determined by the radiation density, can better be viewed as rapidly rising to 
V 1//4 = /c 1//2 (iV) ~ 10 8 GeV then slowly falling to the final reheat temperature @ 
during the reheating process. This reheating gives entropy to the Universe. Non- 
perturbative effects can produce particles with masses up to the potential height, 
i.e. m < Vq^ 4 ~ 10 8 GeU (preheating). The preheating and reheating process in 
this model is quite complicated, but the essential physics is as follows. To begin 
with the potentially problematic axions and gravitinos are not over produced during 
preheating [IT], |2£|. Higgs scalars are copiously produced through preheating via 
the couplings A and k to the oscillating inflaton fields. Although the neutralinos are 
produced in Higgs decays via preheating, once the Higgses decay they go into thermal 
equilibrium, and subsequently freeze out while the universe is radiation dominated, 
similar to the usual hot big bang scenario. However, for a range of parameters the 
singlino is lighter than the lightest neutralino, and in this case after freeze-out the 
lightest neutralino decays into the singlino thereby reducing the LSP relic density by 
the ratio of their masses. 

Recently it has been realised that reheating in all hybrid models, including the 
SUSY motivated ones of interest to us here, goes through very effective tachyonic 



preheating fL6 |. As a result the stage of the background oscillations of the scalars 



around the minimum of the potential will never be reached. This picture is dra- 



matically different from the early papers on the reheating in hybrid models [14, 15 
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and will probably affect the results in [IT]. On the other hand the new picture of 
preheating implies that exciting one field (for example or N) is sufficient to rapidly 
drag all other light fields with which it interacts into a similarly excited state. This 
strengthens our claim that the Higgs doublets which interact with iV will be efficiently 
preheated. 

Once the Yukawa couplings are included in the superpotential, right-handed sneu- 
trinos are also expected to be produced during the initial period of preheating via their 
couplings to the Higgs doublets, and decay out-of-equilibrium into leptons and Hig- 
gses giving rise to leptogenesis. We have already remarked that, unlike the usual hot 
big bang scenario, the out-of-equilibrium condition is automatically satisfied during 
reheating, and furthermore the production mechanism of right-handed neutrinos is to- 
tally different. In the standard scenario the baryon asymmetry is given by Yj, ~ dt\j g* 
where e\ is the lepton number asymmetry produced in the decay of the lightest right- 
handed neutrino of mass Mi, g* counts the effective number of degrees of freedom 
(for the SM g* = 106.75 while for the Supersymmetric SM g* = 228.75) and d is the 
dilution factor which takes into account suppressions from either the couplings being 
too small to thermally produce right-handed neutrinos, or too large to satisfy the 
out-of-equilibrium condition. Typically d <C 1 except for a narrow range of couplings 
| |20| | . In the inflationary picture of reheating outlined above, the baryon asymmetry 
is given by Yj, ~ 76i(cVo) 1//4 / 'M\ where c is the fraction of the total vacuum energy 
converted into right-handed neutrinos due to preheating, and 7 accounts for dilution 
due to entropy production during reheating. Since E\ ~ 10 -6 (Mi/10 10 (jeV') |21] we 
find that Y b ~ 10~ 8 7(c) 1//4 , apparently independently of Mi (although c will depend 
on Mi, for instance c = for Mi > 10 8 GeV.) 

Solving the Boltzmann equations for a particular choice of parameters in this 
model the densities of the neutralinos, radiation, relativistic axions and baryons were 
calculated at reheating time, defined as the time at which the oscillating singlet energy 
density rapidly decayed to zero J7|. This time represents the start of the hot big bang. 
The important point to emphasise is that at this time tRH, for a given model the 
Boltzmann equations enable us to calculate the energy density of the different types 
of matter p ma tter(tRH), as well as the radiation energy density p 7 (t_R//). In the present 
model the details of this calculation are discussed in 0, and a similar calculation 
may be performed for any other model. 
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4 How To Calculate the Size of the Dark Energy 
Density in Supersymmetric Inflation 



We now turn to the question of dark energy, which is not addressed by our super- 
symmetric particle physics based model of inflation. The origin of DE might be a 
traditional cosmological constant with equation of state w = — 1 or some time- varying 
smooth energy (quintessence) with — 1 < w < —0.6 where the upper bound is from 
current observations, and may be in conflict with some quintessence models ||30|| . f\ 



Quintessence models assume a zero cosmological constant and add some arbitrary 
field to account for DE. From the point of view of our inflation model, the simplest 
possibility is to assume that at the global minimum (after inflation) the height of 
the potential is not zero but about 10 _3 eV. This possibility, which just corresponds 
to a standard cosmological constant with w = —1, can be arranged (though not 
explained) by tuning V in Eq.£|; V has to be tuned in any case to give a zero cos- 
mological constant, so this possibility requires no additional tuning to that already 
required in the model. In the no-scale SUGRA model || Vq arises from the mod- 
uli fields in the string theory, and is determined by the non-perturbative physics of 
moduli stabilisation which is not yet understood from a fundamental point of view 
but may nevertheless be parametrised. Of course such a proceedure raises the cosmic 
coincidence question: why should we have p DE ~ p ma tter at the present epoch? Until 
the cosmological constant problem is resolved, there is no way that this question can 
be answered. We reject recent claims to the contrary which are based on setting the 
cosmological constant to zero by hand to start with, since there is always the danger 
in this approach that one has thrown away the baby along with the bathwater. In the 
absence of anything better, some authors have turned to anthropic arguments, but 
many anthropoids reject this approach also as long as alternatives may be possible. 

Is there anything that we can say about DE at the current time from the per- 
spective of our supersymmetric particle physics based model of inflation? Perhaps 
surprisingly the answer is positive: we shall show that we can deduce the present 
day value of dark energy density from the model, together with the measured CMB 
temperature and the Hubble constant, even if the model does not yet specify the 
physics of dark energy! 

A key point of our approach is that a supersymmetric particle physics based model 
of inflation enables us to calculate (in principle at least) the (energy or number) den- 
sities of all forms of radiation and matter (but excluding dark energy) at some early 

6 Note that for a scalar field w = p/p = (KE - PE)/(KE + PE) and if the kinetic energy (KE) 
is small compared to the potential energy (PE) then w is negative, w < — 1 would require negative 
KE which corresponds to the scalar being a ghost field, and a loss of unitarity fl3l|. 
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time tun after inflation and reheating has taken place, corresponding to the start of 
the standard hot big bang. For simplicity we consider only one type of matter energy 
density Pmatteritiw) (which may readily be obtained from the calculated number den- 
sity) and radiation energy density p 7 (£r#). The argument may be straightforwardly 
generalised to the case of several components of radiation and matter. Now, using 
the equations of the standard hot big bang, we wish to obtain their values at the 
present time t , p 7 (to) and Pmatterito)- Without further information this is impossible 
since we need something to tell us when the present time t is, and moreover the 
model does not specify either pDEitRn), or its equation of state, both of which will 
influence the evolution of the universe. Therefore let us input into our analysis the 
present day observed CMB temperature T = 2.725K, which corresponds to a photon 
density p 7 (to) — (2.115 x 10~ 4 eV) 4 , a photon number density n 7 = 410cm~ 3 , and, 
assuming three families of light neutrinos, an entropy density s = 7.04n 7 . Then, 
ignoring additional sources of entropy between t^H and t (such as electron-positron 
annihilation), since we know the equations of state for photons and matter, p 7 ~ R~ 4 
and Pmatter ~ R~ 3 where R is the scale factor of the universe, using the initial values 
of p^rh) and p ma tter(tRH) predicted by the model and the present value of p-y(t ) 
from observation, we find p matt er(to) ~ Pmatter (t RH ) [p 7 (*o)/p 7 3/4 • We emphasise 
that this determination of p m atter{to) is independent of the unknown dark energy. Al- 
lowing for entropy production, which will increase the photon energy density relative 
to the matter energy density, it is usually convenient to consider the ratio n mat t er /s 
which is equal to the number of particles of each species per comoving volume. From 
the model we can calculate n matter /s at and by particle number conservation the 
value of this ratio at the present time t is unchanged. Using the present value of s 
(above) we therefore immediately find n maUer (t ) and from the mass of the particle 
type we readily find p m atter(to), again independently of the dark energy. 

Once we have obtained p m atter{to), from a combination of our model calculation 
and the observed CMB temperature, as outlined above, we now use the observed 
Hubble constant H , or equivalently the present day critical density p cr u, to con- 
vert Pmatter(t ) into the various ^matter = Pmatter (h) I ' Pcrit- Once Vt matter is predicted 
within some supersymmetric particle physics based model of inflation, supplemented 
by measurements of the CMB temperature and the Hubble constant, then it is clear 
that Qde is also predicted to be 

QdE — 1 — ^matter (13) 

Thus a model of inflation that is capable of predicting fl ma tter using measurements 
of the CMB temperature and the Hubble constant, is also capable of predicting Vt DE 
from Eq.|i~3|. This sum rule was written down in ref. [g^], including a curvature term 



and it was discussed there how to determine each of the components VLde and VL matter 
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from observation. What we are saying here is quite different from the empirical 
approach to determining the components of this equation discussed in ref.[|32|], and 



should not be confused with it. To begin with we are assuming inflation, so that 
the curvature contribution is zero. Secondly we are only taking two inputs from 
observation, namely the CMB temperature and the value of the Hubble constant. 
Given these inputs we have shown how an inflation model allows us to then calculate 
flmatter, and hence deduce Qde from Eq.[l3|. 

At first sight our result appears surprising: how can we have deduced Qde from 
a model in which the dark energy is unspecified? In order to answer this it is useful 
to compare two slightly different models of inflation, one in which the dark energy 
density is zero and one in which it is non-zero, but which otherwise predict identical 
values of p-y(t RH ) and PmatteritRn) (in our example this just corresponds to tuning V 
slightly differently in the two cases leaving all the other parameters unchanged.) In 
both cases this will result in the same values of pmatteriio), once p 7 (to) is inputted. 
The only difference is that the Friedmann equations, with zero curvature term due to 
the inflation assumption, will determine two different values of the Hubble constant, 
corresponding to two different values of critical density. The first is given by p cr it\ — 
Pmatter{to), and the second involving the sum of two contributions p cr it2 = Pmatter(to) + 
Pde{^o). Since we input the Hubble constant from observation, we know the true 
value of p cr u in our universe, and so we can discriminate between the two cases from 
a measurement of the Hubble constant. More generally, it is clear that, once the 
correct supersymmetric particle inflation model is known, and the present day value 
of pmatterito) is calculated from it (using the CMB temperature as input), that the 
present day value of Pde(^o) may be deduced from the Hubble constant H which 
is telling us information about dark energy by telling us the critical density. From 
this example it is clear that our argument gives us no new insight into the cosmic 
coincidence question, since a universe without dark energy would simply correspond 
to having a different Hubble constant. 

Once the importance of the Hubble constant Hq in our argument is realised, the 
next question is whether our argument contains any content at all? The answer must 
be yes, since the conclusion relies on non-trivial information coming from the model, 
namely the initial condition for p ma tter(tRH) and p 7 (t,R_ff) without which it would be 
impossible to find p m atteriio) from the CMB temperature alone, and without p ma tter{to) 
it would be impossible to deduce Pde(^o) from the Hubble constant H . A related 
question, is whether our argument is actually useful in practice, given that at the 
present time we do not know the correct model, and even if the model were known 
and all the parameters of that model were accurately specified, we still would need 
to know the physics of preheating and reheating very well. Also when the argument 
is generalised to all the forms of matter and radiation we would need to have a 
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good understanding of the physics of baryogenesis and a way of calculating the axion 
misalignment angle, and so on in order to be able to specify the present day relic 
densities of all the component forms of matter. One could criticise the argument 
on the grounds that the accuracy of the deduced density of dark energy is therefore 
limited by the accuracy with which the matter density can be calculated. While 
this is true, it would seem remarkable to us to suppose otherwise: while it would be 
nice to be able to calculate the DE density much more accurately than the matter 
density, this possibility hardly seems very likely. What our argument gives is a way 
of calculating Qde to the same precision as flmatter, and this we believe is the best 
that one can hope for. 

Whether the DE is due to scalar fields, or an incomplete vanishing of a cosmo- 
logical constant (corresponding perhaps to the universe ending up temporarily in a 
vacuum state slightly above a global minimum at zero), perhaps some of the param- 
eters that determine it will be related to parameters that also determine the forms of 
matter. In the present paper we do not specify any particular model for DE and so 
we must therefore rely on observation to determine the equation of state for the DE. 
If observation eventually tells us that w — — 1 and that the DE is equivalent to a cos- 
mo logical constant, then it will be a tremendous challenge to theorists to explain this 
(see for example the approach of Bastero-Gil, Mersini, and Kanti |33|| ). Explaining a 
vanishing cosmological constant is already proving a very difficult question for string 
theorists, and explaining a very small one does not apparently make this any easier. 
In this case it is possible that the DE question will be around for a long time. In the 
meantime progress may be forthcoming on determining the supersymmetric particle 
inflation model and in determining its parameters and in being able to use those 
parameters to calculate the relic matter densities with increasing accuracy. In such a 
scenario, the one consolation may be that our argument enables one to then calculate 
the size of the DE density (i.e. the cosmological constant) even in the absence of any 
theory of it. 



5 Summary and Conclusion 

We have considered the class of supersymmetric inflation models in which inflation 
occurs at the intermediate PQ symmetry breaking scale. Such models are better 
motivated than GUT scale inflation models which face the problems with magnetic 
monopoles and domain walls. In intermediate scale supersymmetric inflation the 
same theory which is responsible for inflation is also responsible for the solution to 
the n problem, and the strong CP problem. As a consequence one would generally 
expect CDM to contain an axion component in addition to an LSP component. The 
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LSP itself need not be the lightest neutralino, but may be a singlino associated with 
the singlet fields which control hybrid inflation and resolve the [i problem. The 
present day relic densities of the CDM components comprising LSPs and axions may 
be calculated in a given model, using the observed CMB temperature and Hubble 
constant, although the axion density will be subject to the usual uncertainties of the 
unknown misalignment angle. 

Once right-handed neutrinos are added, as the recent confirmations of neutrino 
masses suggests that they should be, then the possibility of baryogenesis via leptoge- 
nesis seems well motivated, and then the baryon density may be calculated in a given 
model. The various relic densities are in principle calculable in a given model, and 
are related to each other and to other phenomena since the number of parameters 
involved is generally smaller than would be the case without a theory. The same pa- 
rameters control cosmology on the one hand and collider physics on the other hand. 
For example common soft SUSY breaking parameters are involved in both inflation 
and collider physics. In order to illustrate all these ideas we have described an ex- 
plicit intermediate scale supersymmetric inflation model which already exists and is 
quite well studied in the literature || [7|, and many of the general ideas above are 
made very explicit in the model, for example the role of the underlying parameters 
in determining different phenomena is demonstrated for this model in Table 1. 

Over the next few years there will be considerable progress in cosmology from the 
Map and Planck explorer satellites, and in SUSY from the Tevatron and LHC We 
believe the time is ripe for a new closer synthesis of SUSY and inflation, and that 
the most promising scenario will involve these theories meeting at the intermediate 
scale. We have shown that in this case one may hope to relate different phenomena 
in cosmology and in particle physics in a much closer and more predictive way than 
ever before. Finally we have made the original observation that, given the value of 
the CMB temperature and Hubble constant from observation, an intermediate scale 
supersymmetric inflation model allows the present day matter relic density to be 
calculated, and hence the present day DE relic density to be determined from Eq.p~3| 
even in the absence of any theory of DE. 
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